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DETERMINATION OF MEAN ORBITAL ELEMENTS FROM 
INITIAL CONDITIONS FOR A VINT! BALLISTIC TRAJECTORY 


by 

Harvey Walden 
Goddard Space Flight Center 


INTRODUCTION 

The determination of a set of constants of motion for an orbital satellite theory is a commonly 
encountered problem inasmuch as mathematical theories of satellite motion are generally given as 
functions of mean orbital elements rather than osculating elements. In practice, it is the set of initial 
conditions of position and velocity components at a given epoch time that is readily available, such as 
those arising from nominal conditions for an orbit insertion maneuver or as the output of a stepwise 
numerical integration technique for trajectory prediction. Such initial conditions are readily con- 
verted, by means of the Keplerian two-body transformations, into osculating orbital elements, but the 
problem of producing mean elements for use as constants of motion in an analytic development re- 
mains. 

A method has been proposed (Reference 1 ) for the determination of mean orbital elements di- 
rectly from initial conditions or Keplerian osculating elements for the spheroidal theory of satellite 
orbits developed by Vinti. The spheroidal theory (References 2 and 3) provides an algorithm for the 
calculation of an accurate reference orbit for any drag-free satellite that moves in the gravitational field 
of an axially symmetric oblate planet. As applied to the actual gravitational potential of the earth, the 
reference orbit accounts exactly for the effects of all zonal harmonic terms in the series expansion of 
the geopotential through the third term, and it accounts for the major portion of the fourth zonal har- 
monic as well. The spheroidal theory is applicable to all bounded orbits of arbitrary inclination and 
eccentricity. The method for the determination of mean orbital elements for Vinti’s satellite theory 
was applied (Reference 1 ) to an actual trajectory corresponding to a near-earth satellite orbit of me- 
dium inclination and moderately high eccentricity that remained above the portion of the atmosphere 
that induces appreciable drag effects. 

A recent paper (Reference 4) describes a method for the utilization of differential coefficients to 
fit orbital parameters to assumed initial position and velocity vectors that represent a ballistic trajec- 
tory. This method adopts an earlier spheroidal satellite theory by Vinti (Reference 5) which does not 
have the advantage of incorporating the effects of the third zonal harmonic term into the spheroidal 
potential. Also, the inversion of the integrals of motion in the form suggested by Izsak (Reference 6) 
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for the solution of Vinti’s dynamical problem is utilized for the coordinates, and the equations for the 
velocities are those given by Borchers (Reference 7). However, the method of determination of the 
Izsak-Borchers orbital elements from initial conditions is nearly identical to the method proposed ear- 
lier (Reference 1) in that both methods are iterative procedures involving a first-order Taylor’s series 
expansion at epoch time that is dependent upon partial derivatives that assume the form of differential 
coefficients. The minor differences between the two iterative procedures shall be discussed further 
below, but it is to be emphasized that the methods are substantially the same. 

The purpose of the present paper, a sequel to the earlier results (Reference 1), is to demonstrate 
the feasibility of applying the method as originally presented to ballistic trajectories. Such trajectories 
are elliptic (or circular) segments of orbital arcs which eventually intersect the earth’s surface (i.e., 
their “perigee” heights are less than unity when measured in units of earth radii). These trajectories 
may correspond to the free-flight portion of an ascent to satellite orbit or to a spacecraft reentry into 
the atmosphere. 

DETERMINATION OF MEAN ORBITAL ELEMENTS 

The constants of the motion q i (i= 1,2, . . . , 6), which are the mean orbital elements for Vinti’s 
spheroidal theory of satellite motion, include the following: the semi-major axis a\ the eccentricity e\ 
a parameter S corresponding to sin 2 i (where i is the inclination of the orbital plane to the Equator) in 
Keplerian two-body motion; and three parameters jSj, (3 2 , and /S 3 , which correspond to the negative of 
the time of passage through perigee t, to the argument of perigee w, and to the right ascension of the 
ascending node £2, respectively, in the reduction to Keplerian motion. The given initial conditions for 
a satellite orbit are generally provided in the form of rectangular inertial position components x 0 , y Q , 
and z Q and velocity components * 0 , y Q , and z Q , specified for a particular epoch time t Q . For the case 
in which the given initial conditions are provided in the form of Keplerian osculating orbital elements, 
namely, a, e, i, r, co, and £2, transformation to inertial Cartesian position and velocity vector compo- 
nents may be readily achieved by means of the usual two-body transformation equations. 

The problem of the determination of the proper set of mean orbital elements from the given ini- 
tial conditions has been approached by Vinti (Reference 8) through a method of factorization of two 
quartic polynomials that arise in the inversion of the integrals of motion. This factorization is carried 
out iteratively, beginning with a zeroth-order solution that corresponds to Keplerian two-body motion. 
A set of four nonlinear equations that results is solved by a method of successive approximations car- 
ried through second order in the oblateness parameter J 2 lp 2 , where J 2 is the coefficient of the second 
harmonic term in the infinite series expansion for the geopotential and p is the semi-latus rectum of 
the orbit measured in units of earth radii. The second-order transformation equations for the orbital 
elements used in the equations of the final solution for the satellite coordinates and velocities (Refer- 
ences 9 and 10) are provided explicitly in Reference 8. To obtain the solution to the nonlinear system 
to arbitrarily high order, a method that applies the Newton-Raphson iteration scheme has been pro- 
posed by Borchers (Reference 7). 

An alternative method for the determination of mean orbital elements for Vinti’s spheroidal satel- 
lite theory that is based upon differential corrections applied to position and velocity residuals at 
epoch time has already been specified (Reference 1). A summary of this method appears to be in 
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order here. This method is capable of determining mean elements directly from initial conditions, 
eliminating the need for numerical factorization through successive approximation, and it has no con- 
nection to the traditional differential correction of satellite orbits that utilizes observational data. 

If the Keplerian osculating orbital elements a, e, i, t, go, and £2 are adopted as the constants of the 
motion in Vinti’s spheroidal satellite theory, then the rectangular inertial position and velocity vectors 
predicted analytically by the theory for the epoch time t Q may be denoted x, y, and z and x, y, and z. 
These components will differ from the initial conditions because the spheroidal theory incorporates an 
earth model that is considerably more sophisticated than the Keplerian elliptic model. If one assumes 
that the required corrections to the orbital elements are sufficiently small so that their squares and 
higher powers may be neglected (as is traditionally the case with linear approximations), the residual 
differences in coordinates and velocities may be expressed by a truncated Taylor’s series expansion: 



Here the predicted values x, y, and z and x, y, and z are considered functions of the six independent 
variables q. (i = 1,2,..., 6) that are to be improved by the respective additive increments Aq j (i = 1 , 

2, . . . , 6). Note that the time t 0 has been omitted as an independent variable inasmuch as it remains 
constant throughout. The 36 differential coefficients appearing in the Taylor’s expansions have been 
evaluated explicitly (References 1 and 1 1) and retain the accuracies of the original Vinti orbital theory. 
The above simultaneous system of six linear algebraic equations admits a solution for the six unknowns 
A q. that are used to correct the orbital elements: 

qi^q.+Aq., 1 * 1 , 2 ,..., 6 . 

The corrected (primed) orbital elements are then used in the spheroidal theory to predict analytically 
a new set of position and velocity vectors at epoch time and to reevaluate the differential coefficients. 
The process of improving the elements is then continued iteratively until the absolute values of the 
position and velocity residuals at epoch time, U 0 - x|, |y 0 -y\, . . . , |z 0 - z|, reach some sufficiently 
small predetermined values. 

When the method outlined above is compared with the method presented by Allen (Reference 4), 
striking similarities can be seen. A truncated Taylor’s series expansion is given by Equations 54 and 
55 of Reference 4 which, after elimination of the apparent typographical errors, is 
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where 


« dB. 

^^ + T,w Ay ‘' 

i= 1 ' 

A%= 7,(13,.) -7,0^). 

In this Taylor’s expansion, j3 ; . = ^.( 7 .) is a component of a matrix consisting of the six Izsak-Borchers 
orbital elements. These are analogous to the mean orbital elements used in Vinti’s spheroidal theory, 
but without the inclusion of the third zonal harmonic in the reference geopotential (References 5 and 
8 ). The parameters 0? ate the Keplerian two-body orbital elements that one obtains by fitting to the 
initial vectors, having set the principal oblateness parameter^ equal to zero in the Izsak-Borchers 
equations of motion. Also, y. represents a component of a matrix which consists of the six oblate 
spheroidal coordinates and their time derivatives that are used in the solution of the Vinti dynam- 
ical problem. The differential coefficients d^./dy j appearing in this form of the Taylor’s series are in 
some sense “inverses” of the previous differential coefficients, inasmuch as they consist of partial de- 
rivatives of the orbital elements with respect to the coordinates and velocities rather than the reverse. 
However, the differential coefficients are evaluated by Allen for the simplified model of a spherical 
earth only and do not include any oblateness effects. In Allen’s version of iterative improvement, the 
mean orbital elements are improved indirectly through the residuals Ay t in the oblate spheroidal coor- 
dinates and velocities, i.e., the differences between the “exact” values 7 f (/J.), based upon the latest cor- 
rected elements, and the approximate values 7 ( (j3?), based upon the elements of the previous iteration. 
Note also that the six Taylor’s series expansions may be solved individually as independent equations 
rather than by Gaussian elimination as a simultaneous system of six linear algebraic equations. 

NUMERICAL APPLICATIONS TO A BALLISTIC TRAJECTORY 

The iterative method of determining a set of mean orbital elements from initial conditions for 
Vinti’s spheroidal satellite theory has already been applied (Reference 1) to an actual satellite orbital 
trajectory that had an inclination of 46 degrees, an eccentricity of 0.24, a period of 195 minutes, and 
apogee and perigee heights of 4600 and 1300 statute miles, respectively. The same method is now to 
be applied to a ballistic trajectory represented by the initial conditions given in the first column of 
Table 1 as inertial rectangular coordinates and their time derivatives. These data are precisely those 
specified by Allen (Reference 4) with the velocity components converted from units of earth radii per 
second, as presumably adopted by Allen, to earth radii per canonical time unit— the latter defined to 
be 806.823 seconds, as indicated in the footnote to Table 1. The second column of Table 1 provides 
the osculating Keplerian orbital elements that correspond to the initial conditions. These classical 
Keplerian elements were obtained through use of the two-body equations of motion, with only the 
central term in the geopotential considered. The third column of Table 1 displays the converged val- 
ues for the mean set of orbital elements for Vinti’s satellite theory obtained after two iterations. A 
measure of the degree of improvement in the orbital elements is provided by Table 2, which displays 
the residuals in position and velocity components at the epoch time. These residuals were obtained by 
the use of the osculating Keplerian elements initially and then by the use of the mean Vinti elements 
obtained upon convergence after two iterations. As an indication of the convergence speed of the itera- 
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tive method, Table 3 presents the residuals in position and velocity at each iteration, where iteration 0 
corresponds to the use of Keplerian elements. 

The results of an improvement in the mean orbital elements for Vinti’s satellite theory, following 
the use of iterative factorization of the quartics through second order, are summarized in Table 4. The 
initial conditions for the ballistic trajectory are identical to those included in Table 1 . The second 


Table 1 —Determination of mean orbital elements from initial conditions. 


Initial Conditions* 

Osculating 
Keplerian Elements* 

Mean 

Vinti Elements* 

x = 0.4650 

a = 0.78809935 

a = 0.78675654 

y = 0.7254 

e = 0.59910027 

e = 0.60188516 

z = 0.6010 

sin 2 i = 0.37911289 

5= 0.37863865 

x = 0.7915 

-t = 1.1628914 

= 1.1638504 

y = 0.03026 

CO = -0.66925097 

0 2 = -0.69029205 

t = 0.08673 

12 = 3.0392151 

j3 3 = 3.0260546 


*The position components x,y , and z are in earth equatorial radii (e.r.), and the velocity components x,y, and z are in e.r./c.u.t, where 
1 canonical unit of time (c.u.t.) is equal to 806.823 seconds. The semi-major axis a is in e.r.; the time of perigee passage r and are in 
c.u.t.; and the argument of perigee cd, the right ascension of the ascending node 12, P 2 , and are in radians. 


Table 2— Magnitude of residuals in position and velocity components. 


Residual* 

Initially 

(osculating Keplerian elements) 

Upon Convergence 
(ipean Vinti elements) 

|Ax| 

25032.81 

0.33 

|Ay| 

9598.04 

0.14 

|Az| 

27478.12 

0.24 

|Ax| 

3282.7721 

0.0236 

|Ay| 

5.2298 

0.0328 

|Az| 

6006.7360 

0.0044 


*The position residuals are in meters, and the velocity residuals in centimeters per second. 


Table 3— Convergence of the computed position and velocity to the initial conditions. 


Iteration 

Position Residual* 

Velocity Residual* 

0 

0.006 019 025 

0.008 659 108 

1 

0.000 085 298 

0.000 064 523 

2 

0.000 000 068 

0.000 000 051 


*The position residual is defined by y/ (Ax) 2 + (Ay) 2 + (Az)^ and is given in e.r., and the velocity residual is defined by 
y/ (Ax) 2 + (Ay) 2 + (Az) 2 and is given in e.r./c.u.t. 
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column of Table 4 provides the mean orbital elements for Vinti’s satellite theory determined by iter- 
ative factorization of the quartic polynomials in the integrals of motion carried through second order 
in the earth’s oblateness parameter. The third column of Table 4 gives the converged values for the 
improved set of mean orbital elements, again obtained after two iterations. Table 5 shows the residuals 
in position and velocity components obtained by use of the factored elements initially and then by use 
of the mean elements obtained upon convergence after two iterations. The convergence of the residuals 
in position and velocity is indicated in Table 6, where iteration 0 now corresponds to the use of the fac- 
tored Vinti orbital elements. Note that “near-convergence” is attained after only a single iteration. 


Table 4— Improvement of mean orbital elements after use of second-order factorization. 


Initial Conditions* 

Factored 

Mean 

Vinti Elements* 

Vinti Elements* 

x = 0.4650 

a = 0.78675646 

a = 0.78675653 

y = 0.7254 

e = 0.60188502 

e = 0.60188516 

z = 0.6010 

S = 0.37864354 

S = 0.37863868 

x = 0.7915 

= 1.1638503 

0! = 1.1638504 

y = 0.03026 

0 2 =-0.69029421 

0 2 =-0.69029207 

z = 0.08673 

0 3 = 3.0260707 

0 3 = 3.0260546 


*See footnote to Table 1. 


Table 5— Magnitude of residuals in position and velocity components. 


Residual* 

Initially 

(factored Vinti elements) 

Upon Convergence 
(mean Vinti elements) 

|Ax| 

91.43 

0.21 

|Ayl 

34.02 

0.05 

|Az| 

27.70 

0.10 

|Ax| 

0.4830 

0.0059 

|Aj>| 

10.8230 

0.0121 

|Az| 

1.3488 

0.0103 


*See footnote to Table 2. 


Table 6— Convergence of the computed position and velocity to the initial conditions. 


Iteration 

Position Residual* 

Velocity Residual* 

0 

0.000 015 900 

0.000 013 810 

1 

0.000 000 107 

0.000 000 083 

2 

0.000 000 037 

0.000 000 021 


*See footnote to Table 3. 
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II 


The conclusions reached in the earlier study (Reference 1) are also valid in the numerical applica- 
tion of the iterative method to a ballistic trajectory. The converged values of mean orbital elements 
for Vinti’s spheroidal satellite theory, as determined from initial conditions, are virtually independent 
of whether or not the process of second-order factorization of the quartics is utilized prior to the iter- 
ative Taylor’s expansion. This is seen from the nearly identical values for the mean Vinti elements pre- 
sented in Tables 1 and 4. Hence, the iterative method of determining mean orbital elements for a 
Vinti ballistic trajectory may be used as a valid alternative to the factorization procedure. However, if 
second-order factorization is applied, then the mean orbital elements are corrected (through subsequent 
application of the iterative improvement method) only by increments of the third order in the oblate- 
ness parameter. Convergence of the iterative fitting to initial conditions in position and velocity is ex- 
tremely rapid, and decreases in the residual components to very close tolerances are achieved, both 
with and without use of factorization to determine initial elements. 
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Appendix A 


In the interest of thoroughness, efforts were undertaken to duplicate the numerical calculations 
of Allen (Reference 4). The following sequence of tabular displays shows that this effort was, at best, 
only partially successful. In Table A1 appear the adopted ballistic initial conditions and, in the second 
column, the Keplerian orbital elements associated by Allen with the initial conditions. The latter are 
seen to differ substantially from the Keplerian orbital elements given in Table 1 . The third column 
displays the converged values for the mean Vinti elements, obtained after three iterations. 

Table A2 shows the residuals in position and velocity components obtained by use of Allen’s 
given osculating Keplerian elements initially and then by use of the mean Vinti elements obtained 
upon convergence after three iterations. The convergence of the residuals in position and velocity is 
indicated in Table A3, where iteration 0 now corresponds to the use of Allen’s given Keplerian ele- 
ments. Finally, the converged values for the mean Vinti elements, obtained after three iterations, 
shown in the third column of Table A1 may be contrasted with Allen’s converged mean orbital ele- 
ments, obtained after four iterations. After the required manipulations described in the footnote to 
Table A1 are performed, the latter appear as follows: a = 0.78681 e.r., e = 0.56072, S = 0.36196, 

(3j = 1.22172 c.u.t., 0 2 = - 0.71421 radian, and 0 3 = 2.94749 radians. 
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Table A 1 -Determination of mean orbital elements from initial conditions using given osculating 

Keplerian elements. 


Initial Conditions* 

Given Osculating 
Keplerian Elements* 

Mean 

Vinti Elements* 

x = 0.4650 

a = 0.78771 

a = 0.78675653 

y = 0.7254 

e = 0.55812 

e = 0.60188518 

z = 0.6010 

sin 2 i = 0.36201 

S = 0.37863866 

x = 0.7915 

-r= 1.23247 

= 1.1638504 

y = 0.03026 

co = -0.73416 

0 2 = -0.69029211 

z = 0.08673 

n = 2.95757 

0 3 = 3.0260546 


♦The units for all variables are as indicated in the footnote to Table 1. The given osculating Keplerian elements were obtained from Allen’s 
data (Reference 4) by the addition of one anomalistic period to the perigee time and by the subtraction of 27 t and 7 r from the argument 
of perigee and the longitude of the node, respectively, at perigee time. 


Table A2— Magnitude of residuals in position and velocity components. 


Residual* 

Initially 

(given osculating Keplerian elements) 

Upon Convergence 
(mean Vinti elements) 

\Ax\ 

20660.69 

0.17 

\Ay\ 

14224.53 

0.14 

|Az| 

20418.46 

0.10 

|AJc| 

2994.5567 

0.0059 

|Aj| 

48612.7070 

0.0180 

|Az| 

5396.8830 

0.0015 


♦See footnote to Table 2. 


Table A3— Convergence of the computed position and velocity to the initial conditions. 


Iteration 

Position Residual* 

Velocity Residual* 

0 

0.005 071 016 

0.061 987 744 

1 

0.003 066 079 

0.000 820 141 

2 

0.000 010 031 

0.000 005 074 

3 

0.000 000 037 

0.000 000 024 


♦See footnote to Table 3. 
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